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ABSTRACT 

The starburst phenomenon can shape the evolution of the host galaxy and 
the surrounding intergalactic medium. The extent of the evolutionary impact is 
partly determined by the duration of the starburst, which has a direct correlation 
with both the amount of stellar feedback and the development of galactic winds, 
particularly for smaller mass dwarf systems. We measure the duration of star- 
bursts in twenty nearby, ongoing, and "fossil" starbursts in dwarf galaxies based 
on the recent star formation histories derived from resolved stellar population 
data obtained with the Hubble Space Telescope. Contrary to the shorter times of 
3 — 10 Myr often cited, the starburst durations we measure range from 450 — 650 
Myr in fifteen of the dwarf galaxies and up to 1.3 Gyr in four galaxies; these 
longer durations are comparable to or longer than the dynamical timescales for 
each system. The same feedback from massive stars that may quench the flicker- 
ing SF does not disrupt the overall burst event in our sample of galaxies. While 
five galaxies present fossil bursts, fifteen galaxies show ongoing bursts and thus 
the final durations may be longer than we report here for these systems. One 
galaxy shows a burst that has been ongoing for only 20 Myr; we are likely see- 
ing the beginning of a burst event in this system. Using the duration of the 
starbursts, we calculate that the bursts deposited 10 53 9 — 10 57 ' 2 ergs of energy 
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into the interstellar medium through stellar winds and supernovae and produced 
3%— 26% of the host galaxy's mass. 

Subject headings: galaxies: starburst - - galaxies: dwarf — galaxies: evolu- 
tion — galaxies: individual (Antlia dwarf, ESO154-023, UGC 4483, UGC 6456, 
UGC 9128, NGC 625, NGC 784, NGC 1569, NGC 2366, NGC 4068, NGC 4163, 
NGC 4214, NGC 4449, NGC 5253, NGC 6789, NGC 6822, IC 4662, SBS 1415 + 
437, DDO 165, Holmberg II) 



Characteristics of the Starburst Phenomenon 



Starbursts refer to periods of int ense, massive star formation that persist on timescales 
much shorter than the Hubble time (ISearle et al.lll973l ). The massive stars formed during a 
burst can alter the gas dynamics, future star formation, and chemical composition of the host 
galaxy through the stars' ionizing radiation, mass loss, and nucleosynthesis. The evolutionary 
effects of starbursts can be dramatic and may be one of the keys to understanding the possible 
link between the evolution of dwarf irregular galaxies (dlrrs) to dwar f spheroidal galaxies 



(dSphs) (e.g.. lLoose fc Thuanlll986l : iDekel & Silklll986l : ISilk et al.lll987l : iDavies & Phillipps 



19881 ) along with externa l effects such as ram pressure and tidal stripping from intera ctions 
with other systems (e.g., Mayer et al. 2001 al ibi: Ivan Zee et al. 2004 : Mayer et al. 2006 ). The 
strength of these feedback processes and their evolutionary impact depend on the SF activity 
and therefore are largely dictated by both the intensity and the duration of the starburst 
event. 

The duration of starbursts h a ve been previously measur e d in individual galax ies (e.g., 



larger samples of galaxies (e.g., iThornley et al. 



Schaerer. Contini. fc Kunthl 11999k iMas-Hesse fc Kunth Il999l; iTremonti et al.l 120011 ) and in 



2000 



Harris et al.l 120041 ) . The most com- 



monly cited timescale for a starburst is 5 — 10 Myr, comparable to the lifetime of massive 
stars. It is logical to presume from the similarity of these timescales that bursts might 
be "self-quenching". In other words, the massive star formation of a burst produces stel- 
lar winds and supernovae explosions that could disrupt the gas fueling the star forma- 
tion and extinguish the starburst on the timescale of a few to ten Myr. The hypothesis 



1 Based on observations made with the NASA/ESA Hubble Space Telescope, obtained from the Data 
Archive at the Space Telescope Science Institute, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under NASA contract NAS 5-26555. 
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of "sel f-quenching starbursts" has been put forth in a number of papers both observa- 



tional (Schaerer. Contini. fc Kunth Il999 ; Mas-Hesse fc Kunthl 



Tremonti et al. 



1998 



200 



Stinson et al 



Harris et al.ll20041) and theoretical ( ITosi et al 



1999 



Thornlev et al. 



1989 



2000 



Ferguson fc Babul 



20071 ). all of which derive < 10 Myr durations for starbursts. This same 
short timescale has been adopted in a number of models. Ferguson et al. (1998) use a burst 
duration of 10 Myr with a star formation efficiency of 10% in modeling the redshift distri- 
bution of dwarf galaxies. Krueger et al. (1995) assume a duration of 5 Myr when modeling 
the spectr al energy distribut i ons of burstin g dwarf galaxies. The sh ort duration assumed in 
models by iKauffmann et al.l ( 119941 ) and by iHogg &: Phinneyl (119971 ) affects the slope of the 
faint galaxy count in the galaxy luminosity function. 



However, our recent study of three nearby dwarf galaxies (jMcQuinn et al.ll2009f ) showed 
both that starbursts can last much longer (i.e., hundreds of Myr) and that starbursts are a 
global phenomenon that can migrate through a galaxy. The duration of starbursts is not set 
by the lifetime of massive stars in these galaxies. These results are in ag reement with earlier 



suggestions that starbursts can la st much longer than a few to ten Myr (jCalzetti et a 



19971 : 

Greggio et al.lll998l : lMeurerll2000l ). A complimentary study by ISchulte-Ladbeck et al.l ( 20011 ) 
reported that although simulations support fluctuations in SFR on short timescales, long 
periods of quiescence between short bursts produce gaps in the simulated color— magnitude 
diagrams (CMDs) not seen in observational data. 

Starburst driven winds may feed the intergalactic medium (IGM) altering the internal 
kinematics and evolutionary path of the host galaxy. The importance of galactic winds and 



metal— rich gas outflows triggered b y sup ernovae explosions in low mass ga 



axies has been 



much debated in the literature (e.g., 


Dekel & Silk 


1986; 


Skillman et al. 


1989; 


Skillman 


1997; 


Spaans & Norman! 


1997 




Mac Low & Ferrara 


1999; 


Romano et al. 


2006 




Dalcanton 


2007; 


Martin et al. 


2002; 


Tolstov et al. 


2009, and references therein"). If starbursts are sustained 



on timescales of a few 100 Myr for the majority of dwarf starburst galaxies, longer bursts 
may more readily power galactic winds. 

Here, we directly measure the durations of starbursts in twenty nearby bursting dwarf 
galaxies from their star formation histories (SFHs) reconstructed from archival Hubble Space 
Telescope (HST) observations. We present ed the details of the ob servations, data reduction, 
and method of reconstructing the SFHs in lMcQuinn et al.l ( ]2010l . hereafter, Paper I). In this 
paper, we provide a summary of the observations and data processing in £j2j measure the 
duration of the burst in each system in §3j and compare the durations with the dynamical 
properties of each host galaxy in §H The bursts' durations are compared in §5] with the 
5 — 10 Myr timescale often cited. Finally, we calculate the amount of stellar mass and 
energy produced in the burst in $6j Our conclusions are presented in 
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2. The Sample Galaxies and Their Star Formation Histories 
2.1. Observations and Photometry 



Our sample was selected from observations in the HST public archive based on three 
photometric conditions. First, we require both V and / band images of a galaxy in the 
archival observations. Second, we set a minimum photometric depth of the I band obser- 
vations of ~2 mag below the tip of t he red giant branch, which is required t o accurately 
constrain the recent SFH of the galaxy (jDolphinll2002l ; iDohm-Palmer et al.ll2002l ). Third, the 
galaxies are required to lie close enough such that their stellar populations were resolved by 
HST imaging instruments. These properties were chosen to ensure robust reconstruction of 
the SFHs and measurements of the duration. The photometry for one galaxy, SBS 1415+437, 
stopped ~ 0.5 mag short of our depth requirement despite long integration times. We in- 
clude the results on this well— studied galaxy in our sample for comparison purposes. The 
resulting sample of twenty galaxies is listed in Table [TJ 



All photometry used DOLPHOT or HSTphot toolphinl |2000j) on the standard HST 
pipeline processed and cleaned images. The details of the observations and photometric 
processing for eighteen of these ga l axies are presented in Paper I and for two galaxies, 
DDO 165 and Ho II, in lWeisz et al.l (120081 ). As a representative example of the data quality 
and photometry, we present the imaging and corresponding CMD for one galaxy, UGC 9128, 
in the top two panels of Figure [TJ The evolutionary features that drive the reconstruction of 
the SFHs can be readily identified in the CMD including the main sequence (MS) stars, blue 
and red helium burning (BHeB and RHeB respectively) stars, the red giant branch (RGB) 
stars, and the asymptotic giant branch (AGB) stars. The upper MS, the BHeB, and the 
RHeB stars are unambiguous signs of star formation in the past few 100 Myr. The helium 
burning (HeB) stars are of particular interest as a star in this stage of stellar evolution 
occupies a unique space in the CMD; the luminosity of a HeB star has a one— to— one 
correspondence to the mass of the star and, hence, the age of the star. Under the assumption 
of a universal initial mass function (IMF), this luminosity— mass— age relation of the HeB 
stars tightly constrains the SFH at recent times. 



2.2. SFH Method 



We reconstructed t he past rate o f star formation (SFR(t)) for each gal axy by applying a 
CMD fitting technique ( jDolphinll2002l ) with stellar evolutionary models from lMarigo fc Girardi 
( 120071 ) to the stellar photometry. During fitting, the distance and extinction were allowed to 
vary within limited ranges, and the metallicity was constrained to increase with time except 
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in three cases (Antlia, NGC 2366, and Holmberg II). The metallicity in these three systems 
could be fully constrained by the information contained in the deeper photometry reaching 
a magnitude below the red clump. The lifetime SFH results o f the fitting progr am (i.e., 



SFR(t, Z)) are presented in Paper I for eighteen galaxies and in IWeisz et al.l (120081 ) for two 
galaxies (DDO 165 and Ho II). As a representative example of our results, we present the 
complete SFH of UGC 9128 in the bottom left panel of Figure [TJ 

UGC 9128 provides an interesting example of how the features in a CMD translate into 
a SFR (see Figured]). The distribution of stars along the helium burning (HeB) sequences is 
clumpy, suggesting that the SF has been elevated during several distinct episodes interspersed 
with lower levels of SF. We have highlighted stars produced during three periods of elevated 
SF labeled as events (A), (B), and (C) in the CMD. The most recent star formation event (A) 
is traced by BHeB stars, RHeB stars, and MS stars in the CMD. The earlier events (B) and 
(C) are traced by BHeB stars, as well as by RHeB stars that are blending into the red clump 
and MS stars that fall below our photometric depth. These three distinct SF events are 
recovered in the last 500 Myr in the SFH as shown in the final panel of Figure [TJ Event (C) 
contains the lowest luminosity BHeB stars in the three events corresponding to the lowest 
mass (i.e., oldest aged) BHeB stars. This SF event is seen in the SFH as an elevated SFR 
from ~ 250 Myr ago to ~ 400 Myr ago. More recently (i.e., above the luminosity of event 
(C) in the CMD) is a period of relative quiescent SF followed by another increase in SFR 
(the clump of BHeB stars labeled event (B)). This more recent SF event is recovered in the 
SFH beginning ~ 200 Myr ago and lasting ~ 100 Myr. Similarly, event (B) was followed 
by a comparative lull in SF before experiencing another increase in SF from ~ 50 Myr to 
~ 25 Myr ago, event (A). UGC 9128's unique SFH demonstrates how directly the HeB stars 
correlate to SFRs at recent times. 

In Figure [2J we show the SFH over the last 1.5 Gyr for all twenty galaxies. These 
recent SFHs share the common characteristic of elevated SF at recent times, as expected 
for our sample selection. However, there are differences within the sample. For exam- 
ple, some galaxies show a sustained, elevated rate of SF over a few hundred Myr (e.g., 
ESO 154-023, NGC 4068, UGC 4483). Others have significant fluctuations in their SFRs 
on shorter timescales (e.g., UGC 9128, NGC 4214). Importantly, five galaxies show "fossil" 
bursts, for which the SFR has returned to a lower level in the most recent time bins (e.g., 
Antlia, UGC 9128, NGC 4163, NGC 625, and NGC 6789). One galaxy, NGC 6822, shows a 
burst just beginning with elevated levels of SF only in the last 50 Myr. 
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2.3. Photometric Depth and Temporal Resolution of SFHs 

As the photometric depth of each galaxy reaches below ~ 2 mag below the TRGB, 
the recent SFHs in our sample have uniform temporal resolution regardless of the exact 
photometric depth. As demonstrated in the case of UGC 9128 in Figured! the HeB stars 
can be unambiguously age-dated from an optical CMD due to the rapid luminosity and 
spectral evolution of the youngest, most massive stars, all of which are well-sampled in 
our CMDs. Thus, the recent SFHs of starburst galaxies are not only robustly derived, but 
temporal resolution ranging from 5t~ 10 Myr to a few 100 Myr, depending on look-back 
time, can also be achieved universally at recent times (t<l Gyr). 

The time resolution of the intermediate and ancient SFH is, in part, defined by the 
depth of photometry. Between l<t<2 Gyr, the SFRs are anchored by stars in the RGB and 
AGB phase of evolution above and just below the TRGB that h ave been shown to correlate 



to SF in this time period (e.g.. lMateolll998t iGallart et al.ll2005l ). At intermediate times, the 



derivation of the SFRs is complicated by the age- degeneracy of RGB stars. However, the 
SFH can be still constrained by bright AGB stars that lie above the TRGB as these stars do 
not suffer from the same age-degeneracy as the RGB stars, providing the photom etry reaches 



deeper than ~ 2 mag below the TRGB (IGallart et all 120051 : IWeisz et al.1 12008 . see Fig. 4). 



At ancient times, deeper photometry reaching depths of the oldest MS turn-off (oMSTO) is 
needed to break the age- degeneracy of stars in the CMD and hence the co-variance of the 
SFR between time bins. Thus, without deep photometry, larger time bins are justified at 
higher look-back times. 

As the time resolution at older look back times depends on the depth of photometry, we 
measured the inherent time resolution of our data by reconstructing the SFH of synthetic 
photometry with depths corresponding to our deepest (I~ 2 mag), typical (I~ —1 mag) 
and shallowest (I~ —2) data. For each photometric case, we input SF that occurred over 
a short time frame (i.e., log(£t)=0.05) at different look-back times. The results from the 
artificial star tests for our real data were used in the SFH derivation thereby taking into 
account completeness and crowding effects. We present the results of these tests in Figure EJ 
The top panel shows the SF input at 100 Myr, 200 Myr, 450 Myr, and 815 Myr as dashed 
lines for three different photometric depths; the recovered normalized SFRs are shown in 
blue (100 Myr), red (200 Myr), cyan (450 Myr), and magenta (815 Myr). Generally, these 
results show that the SF is well-recovered with the peak SF occurring in the correct time bin 
with some SF spreading to the adjacent time bins. This "spread" defines the inherent time 
resolution of the data. For comparison, we overlaid the time binning we employ in our SFH 
results on these plots as shaded gray regions. The time resolution of the SFH is well-suited 
to measuring how long SF events last at recent times. 



-7- 



3. Durations of Starbursts in Dwarf Galaxies 

3.1. Quantifying Starbursts 

The measurement of the duration of a starburst event requires a definition of "bursting" . 
To establish a quantitative threshold, we apply a metric to evaluate the rece nt SFRs relative 



to the historical SFRs over the past 6 Gyr within a host galaxy, as discussed in M cQ.uinn et al. 
j2009h . We used a birthrate parameter (b= SFR / < SFR >) as defined bv IScalo! jl98fih 



which compares the current SFR with the lifetime average SFR. A normalized SFR removes 
the bias against identifying bursting SFRs in smaller galaxies where the absolute SFR is 
generally lower than in larger galaxies. However, rather than using the lifetime average SFR 
as a baseline, we used the average SFR over the last few Gyr, providing a better comparison 
of the recent SFRs with the more recent evolutionary state of the host galaxy. This choice 
both avoids suppressing b by excluding the higher levels of SF typically expected during the 
initial period of galaxy assembly and leverages the most secure part of the SFH. We use 
our modified birthrate parameter (b recen t =SFR /<S FR> n _6 o vr ) to define a threshold for 



bursting SFRs; a value of b recent > 2 identifies a burst ( IKennicutt et al.ll2005[ ) and b rece nt — 1 
marks the beginning and end of a burst event. In Figure El we show the average SFR from 
the last 6 Gyr (b recent = 1) as a solid red line and twice this average (b recent = 2) as a dashed 
red line. 

The duration of the starburst event can be measured directly as the amount of time 
brecent > 1 once the b recent > 2 threshold has been reached. The time measurement extends 
from the midpoints of the beginning and end time bins. The uncertainties on the durations 
are the width of the longest time bin contributing to the duration. The resulting durations 
of the burst events are shaded in blue in Figure El Fifteen galaxies have starburst durations 
between 450 ± 50 and 640 ± 190 Myr as listed in Table Ho The longest bursts we identify 
are in UGC 9128 and DDO 165, which show bursts lasting 1300 ± 290. The average burst 
duration is ~ 600 Myr. In fifteen cases, the SFR remains above the 6 Gyr historical average 
in the most recent time bin, indicating that the starburst is ongoing. The durations for these 
starburst events are thus lower limits. The average SFRs over the past 6 Gyr, peak SFRs 
and corresponding b recent value during each burst event are also listed in Table [TJ 

Note that the duration measurement is affected by the inherent time binning of the 



2 Note: We measure slightly diff erent durations for NGC 4163, NGC 4068, and IC 4662 than previously 
reported in iMcQuinn et al.l (|2009h . As discussed in Paper I, we constrained the metallicity of the galaxies 
to increase with time, which gives improved constraints on the recent SFH. As a result, the SFH of these 
three systems changed, altering the measured duration of the starbursts by 25% for NGC 4068 and IC 4622 
assuming a double burst, and 50% for NGC 4163. 
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SFH. The temporal resolution over the last 100 Myr is fine enough that we can discern and 
measure fluctuations in the SFR over short timescales and pinpoint the value of b recent every 
20 — 50 Myr or so. In contrast, between 650 — 1000 Myr ago, the inherent time resolution of 
the data reduces to 350 Myr. One can imagine a burst beginning in the latter half of this 
time bin but the b rece nt value is suppressed below unity due to averaging with a lower SFR 
earlier the same time bin. In such a case, this entire time bin would not be considered part 
of the burst event and the measured duration would be underestimated. The reverse is also 
possible. The resolution of time bins has a smaller effect on the final duration measurement 
for bursts contained in the most recent 650 Myr due to the finer temporal resolution. This 
effect is quantified in the uncertainties listed for the durations which increase with larger 
time binning. 

One concern is that our results are in some way dependent on how we bin the old- 
est SF. In order to investigate the effect on our results of dividing the intermediate and 
ancient SF into separate bins, we calculated the duration of the starbursts without distin- 
guishing between intermediate and ancient time bins (i.e., the classical b parameter = SFR 
/ <SFR>jj/ et j me ). Given that Local Group galaxies, whose photometry reaches below the 
red clump, show detect able SF older than 10 Gyr based on red giant population and/or RR 



Lyr stars (lMateolll998l ). we assume the lifetime of the galaxies extend back to our oldest 
time bin. In most cases, including the ancient SF raises the baseline average SFR due to the 
higher rates of SF typical during the initial galaxy assembly period in the early universe. 
Thus, with a higher baseline SFR average, a higher rate of SF is needed to reach the bursting 
threshold of b>2. Nonetheless, only three of the galaxies in our sample (Antlia, NGC 4163, 
NGC 784) did not meet this criterion and would not be considered bursts by this metric. 
The remaining seventeen galaxies are still identified as starbursts. Perhaps surprisingly, the 
duration of the bursts measured with this alternate metric are shortened in only six of the 
seventeen systems, lowering the average duration of the sample from ~600 Myr to ~450 
Myr. Thus, our main scientific conclusion that starbursts have longer durations is robust 
and the uncertain splitting of the time bins is of negligible import. 



3.2. Cases of Complex SFHs 

Our criterion for defining a starburst serve as a prescriptive guide for measuring a burst 
duration but we stress that the criterion is a guide and the distinction between "regular" and 
"bursting" SF is subjective in some cases. In thirteen of the galaxies (e.g., ESO 154 — 023, 
NGC 784, NGC 5253, among others), the SFR profile increases substantially from one time 
bin to another and the duration measurement is easily made from analyzing the b re cent 
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parameter. In NGC 6822, the burst began ~ 50 Myr ago. We measure the lower limit of the 
duration in this system to be 20 ± 13 Myr. Given the much longer durations of the other 
galaxies, this burst is likely just beginning. In six systems the SFH is more complex and 
thus we discuss each system in turn. 

For NGC 2366 and Ho II, the burst is preceded by time bins whose b rece nt is equivalent 
to unity within the measured uncertainties (we refer the reader again to the blue shaded 
regions and red lines in Figure [2]). By our definition, these time bins should be included in 
the duration measurement. Yet, these SFRs continue for hundreds of Myr at roughly the 
historical average SFR. The overall SFR profile suggests that the bursts had not yet begun 
in these time bins. The lower temporal resolution at these older times also affect the b re cent 
calculation. We chose to exclude these time bins from the duration measurement. 

In the third case, NGC 6789, the b recent index drops just below unity approximately 350 
Myr ago after which the SFR increases again pushing the b reC ent close to 2. There is some 
co- variance between adjacent time bins in the SFHs and this short dip in SFR surrounded by 
higher levels of SF suggests that the burst is continuing. The SFR profile indicates that the 
burst event, while weaker than others in our sample, has continued to elevate the SF above 
the average until our most recent time bins. This conclusion is supported by looking at the 
overall SFR profile over the past 1.5 Gyr. The period preceding the burst is a quiescent 
period experiencing very little SF. We chose to include the SF ~ 350 Myr ago as part of the 
burst event which extends the burst duration measurements from 600 Myr ago until the last 
~ 100 Myr. The profile is similar to UGC 6456 where the burst event is clearly defined and 
extended. 

In the fourth case, UGC 4483, the SFR profile experiences a dip ~ 450 Myr ago and 
the b re cent index is equivalent to unity within the measured uncertainties. The time bins on 
either side are elevated reaching b recent values of 13 and 3.9. This suggests that the burst 
has continued although the system experienced a period of lower SF during the larger burst 
event. Likewise, in SBS 1415+437, the SFR dips to zero during one time bin flanked on 
one side with b re cent values of 5.7 and 6.3. We chose to define the burst event inclusive 
of this period of quiescence based on the overall SFR profile. Note that the photometry 
of SBS 1415 + 437 falls ~ 0.5 mag short of our photometric depth requirement, thus the 
ancient SFH is not well-constrained. However, the recent SFH is well-constrained and the 
SFR profile in the last 500 Myr clearly defines a burst. We measure the duration of the 
starburst in SBS 1415 + 437 despite the uncertainties in the past SFR average and include 
it for comparison with others in our sample. 



Finally, the well— studied starburst galaxy, NGC 625 (jCannon et al.l 120031 ). does not fit 
our metric for identifying a starburst. The SFH of NGC 625 is similar to that of Leo A in 
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that the initial SF is suppressed ( ICole et al.ll2007l ) until ~ 3 Gyr ago. The photometry used 
to derive the SFH of NGC 625 reaches nearly three magnitudes below the TRGB giving 
confidence to the low SFRs derived for the 6—14 Gyr time bin. The result of the later 
onset of SF is that the average SFR over the past 6 Gyr is high and the elevated levels 
of SF in the last 500 Myr fall below our starburst threshold. Although the unique SFH 
of NGC 625 does not fit our metric, evidence at other wavelengths classify this system as 
a starburst reinforcing the premise that the identification of starburst (i.e., the difference 
between "regular" and "bursting" SF) is not always clear. For example, if we apply a classical 
birthrate parameter utilizing the SFR average over the lifetime of the galaxy, the recent SFRs 
do meet the starburst threshold. Using the elevated SFRs in the recent history of NGC 625 
to guide the measurement, the duration of this likely starburst event is 450 ± 50 Myr as 
seen in Figure [U much like the majority of others measured in this study. We include it in 
our duration measurements for comparison. Excluding this one unusual system would not 
change our scientific conclusions. 



4. Bursts: SF on Timescales > Few 100 Myr 



The durations of the burst events fall unambiguously on the scale of hundreds of Myr. 
The bursts are sufficiently robust that they are not shut off by an individual flicker of 
SF. Moreover, the longer durations suggest that a starburst event is a galaxy-wide phe- 
nomenon, which likely includes SF in both clusters and in the field regions of the galaxies 
(jMcQuinn et al.ll2009l ). Paper III will specifically investigate the spatial distribution and 
possible migration of SF in the stellar clusters and diffuse regions within the starbursts to 
better characterize the nature of starbursts. 

The burst durations can be compared with the characteristic timescale to communicate 
an event or disturbance across a galaxy. This dynamical timescale can be esti mated as the ro- 



tatio n period for dwarf galaxies as these systems exhibit solid— body rotation (jSkillman et al. 
19881 ) (i.e., tdyn — 2nRe/V where R e is the effective radius of the galaxy and V is the ro- 
tational velocity). We have calculated the dynamical timescale for eighteen^ of the galaxies 
using the major axis radii of the galaxies and the inclination corrected maximum rotation 
curve velo cities. The radii an d velocity measurements were obtained from the HyperLeda 
database (IPaturel et al.l 120031) with the ex c eption of NGC 6822 whose maximum rotational 



velocity was obtained from lWeldrake et al.l (120031 ). The dynamical times calculated are listed 
in Table [1] and are compared with the starburst duration we measured for each galaxy. The 



3 Rotational velocity data were not available for SBS 1415+437 and NGC 6789. 
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ratio of the duration to dynamical timescales for fifteen starbursts range from 0.70 — 4.6 
with a mean ratio of 2.6. Galaxy UGC 9128, with a duration of over 1 Gyr, has a larger 
duration/dynamical timescale ratio of 10; NGC 6822 has ratio of 0.14 as the burst is just 
beginning in this system. The dynamical timescales are plotted as green horizontal bars 
in the SFHs in Figure [2] and can be compared to individual galaxy duration measurements. 
Previous authors have proposed that if starbursts persist only on the characteristic timescale 
of a system or short er, then they w ould likely be self-extinguishing events destroyed by their 



own energy output (IMeurerl 120001 ) . Durations comparable to or longer than the dynamical 
timescales support our conclusion that these events are not self-quenching on the timescale 
of massive star lives as discussed in Paper I. 

The starburst durations we measure are clustered around 450 Myr, as shown in Figure SJ 
It is important to note that 75% of these systems have not yet finished; we are measuring 
lower limits for the burst durations in these systems. The peak of the distribution in Figure H] 
would likely shift to longer timescales if we had a sample of exclusively fossil bursts. The 
distribution of the dynamical timescale is plotted as a dashed red histogram in Figure HJ 
This distribution is similar to the distribution of the burst durations but is shifted towards 
shorter times suggesting the burst process is global in nature. Rapid self-quenching of SF 
would limit starburst durations to ~ 10 Myr, indicated by the dashed blue line. 



5. "Flickering": SF on ~10 Myr Timescales 

Burst durations on the scale of hundreds of Myr differ substantially from the often 
reported timescale of 5 — 10 Myr (see references in JQ). The nearly two orders of magnitude 
difference between these two scales is the result of measuring two different scales of SF 
events; i.e., whether one is measuring the duration of a complete starburst event in a galaxy 
or measuring the duration of a pocket of enhanced SF made up of one or a few stellar 
clusters or associations within a galaxy. Both scales of SF are evident in our SFHs. We 
have highlighted the burst events in blue in Figure |2J Now we draw the reader's attention 
to the most recent time bins with the finest resolution (t ~ 4 — 10 Myr, t ~ 10 — 25 Myr) 
presented in Figure [2j Over these short timescales, there are rapid variations in the SFRs 
seen in many of the galaxies corresponding to local "flickering" of SF on small scales. This 
flickering occurs within a larger envelope of enhanced SF on timescales greater than a few 
100 Myr. The larger starbursts events are likely comprised of numerous individual pockets 
of flickering SF in stellar clusters or associations as well as SF in diffuse field regions and the 
starburst duration is the summation of all these smaller scale events. 

Obviously, small scale increases in the SFR can occur outside of a burst event. In 
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other words, flickering SF and the formation of stellar clusters and associations do not occur 
solely in starbursts. For example, consider the SFH of the Antlia dwarf galaxy in Figure |2j 
The burst event in Antlia ended nearly 200 Myr ago. Yet, we see elevated SF during one 
additional time period (t ~ 10 — 25 Myr ago) where the b rece nt index exceeds 2. This increase 
in SF is temporally isolated from the larger burst event and is likely caused by SF in an 
individual stellar association or small cluster. We do not include this flickering SF as part of 
the burst event as the surrounding time bins have lower levels of SF. One of the implications 
of longer burst durations is that the burst is likely occurring on large spatial scales within 
a galaxy and small increases in SFR are qualitatively different than bursting SF. We will 
explore this in detail in Paper III with a spatial analysis of the SF. Note, we are only sensitive 
to this episode in Antlia because of the fine time resolution in the SFH 25 Myr ago. If this 
isolated flickering event had occurred 200 Myr ago, it would have blended into the SFR 
derived for the larger time bin. 

The variations of SFR on short timescales (5t ~ 5 Myr) are also seen when the SFR 
results in our most recent time bin of 4 — 10 Myr are compared with another starburst 
indicator, Ha emission, as discussed in Paper I. The timescales for our most recent time bin 
and Ha emission overlap but do not completely coincide. As one would expect, the results 
of these two indicators agree in most galaxies in the sample (i.e., where there is significant 
Ha emission, we find b recen t values above two in the 4 — 10 Myr time bin and vice versa). 
Yet in four galaxies, the two indicators disagree suggesting that the SF activity has changed 
rapidly on timescales of a few Myr. This short timescale is the same as what is expected for 
the evolution of massive stars. 

The characteristic timescale to communicate an event or disturbance for "flickering" SF 
is the crossing time for a typical star cluster (t cross = 2R e /a where a is the velocity dispersion 
of the stars). Assuming a typic al radius of a star cluster to be a few pc w ith a velocity 



dispersion of ~ 10 km s 1 (e.g., ISmith fc Gallagher! |2000| ; lLarsen et al.ll2008l ). the crossing 



time is ~ 1 Myr. T his value agrees with the typical dynamical star cluster timescale of a 



few Myr reported by iHeckmanl (120051 ) using a timescales defined by td yn — (Gp)^ 1 ^ 2 . These 
cluster crossing times support the conclusion that the short durations reported for bursts 
are measuring the SF phenomenon within star clusters, associations, or some equivalently 
smaller region of SF. 

The amplitude of the flickering can be used to estimate the stellar mass created in these 
smaller scale SF events. Using a range of SFRs from 0.001 — 0.1 M Q yr _1 typical of the 
sample dwarf galaxies and a timescale of 10 Myr, the total stellar mass produced ranges 
from 10 4 — 10 6 M , which are typical masses of stellar associations and clusters. These 
numbers can be compared to the larger values of 10 7 — 10 9 M Q of stellar mass created during 
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the bursts events of these same galaxies. 



The differences in timescales, dynamical times, and amount of stellar mass produced 
highlight the different physical processes governing local, small scale SF and large scale 
burst events. The flickering SF is governed by local physical conditions (e.g., gas density, 
feedbac k from massive stars, e tc.) and is likely to be temporarily suppressed/quenched 
locally (lOmukai fc Nishil Il999l ; iNishi fc Tashirol 120001 ) . The same feedback from massive 
stars that may quench the flickering SF does not disrupt /quench the overall burst event 
in our sample of galaxies. Thus it appears the non-equilibrium energy output and mass 
tran sfer characteristic of ste l lar winds and supernovae that have been proposed to quench 
SF jFerguson fc Babul! Il998l : iThornlev et all I2OO0I : Istinson et al1l2007l : iKavirai et al.ll2007h 
cannot be assumed universally on galactic scales. 



6. Stellar Mass and Energy Created in a Burst 

Longer duration bursts have a greater impact on the evolution of the host galaxy than 
short bursts. Two interesting parameters to quantify are the percent of stellar mass created 
during a burst and the energy produced by the burst, which can potentially be deposited 
into the ISM. We have calculated the percent of stellar mass created in each burst relative to 
the stellar mass created over the lifetime of each galaxy directly from the SFHs. The results 
range from 3%— 26% of the mass being created during the burst; the individual numbers for 
each galaxy are listed in Table |2j The percentage is an upper limit for SBS1415+437 whose 
ancient SFH remains uncertain due to the limited depth of photometry. The galaxies showing 
bursts with lower absolute SFRs (i.e., Antlia) or SFRs exceeding the past 6 Gyr average in 
only one or two time bins (i.e., UGC 6456, NGC 4163, and NGC 6789) lie at the lower 
end of this range. Galaxies with longer durations and stronger bursts, such as UGC 9128, 
DDO 165, NGC 4214, lie at the upper end of this range. Additionally, we compare the 
total amount of stellar mass calculated form the SFH with the luminosity of the galaxies 
from publis hed absolute magnit udes in the B band, adjusted for extinction from Galactic 



dust maps (ISchlegel et al.lll998l ). The mass derived from the SFH does not incorporate all 
the mass from the optical extent of the galaxies as the field of view in some cases does not 
extend to the outer optical limits. While this technically makes the mass a lower limit, it is 
a small effect as most of the optical area is covered in each system and the amount of mass 
at these outer radii is small compared to the mass in the inner parts of the systems. The 
one exception is NGC 6822 for which the observations do not cover the main extent of the 
stellar populations. The mass— to— light ratios (M*/Lb), in solar units, range from 1.0 ±0.2 
to 8.0 ±1.0. The highest ratio is for Antlia dwarf galaxy, the lowest surface brightness galaxy 
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with an absolute magnitude of —10.14. These M*/Lb ratios are consistent with the expected 
ratios for SF galaxies thus providing an independent check on the derived SFHs. 

We calculated the energy produced in each b urst using our SFHs and STARBURST99 
( jLeitherer et al.lll999t IVazquez fc Leithererl 120051 ) . STARBURST99 is an evolutionary syn- 
thesis code that uses a known SFR and stellar evolution models to predict spectrophotometric 
properties, the energy distribution, and other related properties of galaxies with active star 
formation. Of particular interest here is the ability of STARBURST99 to simulate the energy 
produced by stellar winds and supernovae at a given SFR. Using an instantaneous burst of 
SF with a fixed mass of 10 6 M Q coup led with the same IMF and stellar evolution models 
employed in the SFH code as done in I Weisz et all ( 120091 ). we simulated the energy output 
of SF for low metallicity galaxies with a time resolution of 5 Myr. This energy output for 
a 10 6 M instantaneous burst was scaled by the stellar mass created every 5 Myr for the 
actual observed bursts and integrated over the duration of the burst. The total energy for 
the bursts ranges from 10 53 - 9 — 10 ergs; individual values are listed in the final column 
of Table [2j There is a clear trend that higher luminosity galaxies have a higher amount 
of energy produced in a burst. The significant amount of energy deposited into the ISM, 
particularly in the cases where the energy is > 10 56 ergs, will have an important impact 
on the gas dynamics and future SF activity of the host galaxy. We will study the spatial 
distribution of the SF activity in Paper III and explore whether this energy is deposited in 
concentrated areas or over a significant fraction of the galaxy. 



7. Conclusions 



The starburst events in the the twenty dwarf galaxies studied present bursts lasting on 
the order of hundreds of Myr. The shortest durations measured were ~ 450 Myr, the longest 
were 1.3 Gyr. The majority of the galaxies display bursts with durations ranging from ~ 450 
Myr to ~ 650 Myr but many are lower limits as the SFR is still elevated at the most recent 
times. The inherent time resolution of our data allows us to accurately derive the duration 
of starbursts occurring over the past ~ 1 Gyr. The ratio of the durations to the dynamical 
timescales of the galaxies ranges from 0.70 — 4.6 for sixteen of the galaxies with one outlier, 
UGC 9128, having a ratio of duration to dynamical time of 10. These longer duration times 
rule out the theory that these starbursts are "quenched" by the feedback processes of the 
massive stellar populations created. Shorter durations often cited on the order of 5 — 10 Myr 
are measuring small pockets of enhanced SF or "flickering" that contribute to a burst event 
but do not represent the extant of a burst. The bursts deposit 10 53 9 — 10 57 ' 2 ergs of energy 
into the interstellar medium through stellar winds and supernovae and produce 3%— 26% of 
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the host galaxy's mass. 
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Table 1. SFRs, Burst Durations, and Timescales 



Galaxy 
(1) ' 


<SFR> 6 _q Q yr Peak b rcccnt 
(10~ 3 M yr" 1 ) of Burst 
(2) (3) 


Peak SFR of Burst 
(10- s Mq yr- 1 ) 
(4) 


Duration 
(Myr) 
(5) 


Tdyn 

(Myr) 
(6) 


Dur/T dyn 
(7) 


Antlia Dwarf 


0.15 ± 0.01 


3.5 ±0.9 


0.52 ±0.01 


640±190 


190±20 


3.3±1.0 


UGC 9128 


0.82 ± 0.07 


6.3 ± 1.4 


5.1 ± 1.0 


1300±300 


120±10 


10.±7 


UGC 4483 


0.79 ± 0.06 


14 ±3 


11 ± 2 


> 810±190 


200±23 


4.1±1.0 


NGC 4163 


4.0 ±0.1 


2.9 ±0.6 


12 ± 3 


460±70 


240±20 


1.9±0.1 


UGC 6456 


3.0 ±0.2 


7.6 ±1.1 


23 ±3 


> 570±60 


160±30 


3.5±0.7 


NGC 6789 


3.9 ±0.1 


3.8 ± 1.3 


15 ±5 


480±70 






NGC 4068 


10. ±0.3 


4.7 ±0.3 


46 ±3 


> 450±50 


270±20 


1.7±0.1 


SBS 1415+437 


12 ±2 


12 ±2 


150 ± 10 


> 450±50 






DDO 165 


13 ±1 


6.5 ±0.5 


80. ± 5 


> 1300±300 


710±40 


1.8±0.4 


IC 4662 


9.8 ±0.6 


7.7 ± 1.6 


76 ± 15 


> 450±50 


100 


4.6 


ESO 154-023 


20. ± 1 


6.4 ±0.5 


120 ± 10 


> 450±50 


460 


1.0 


NGC 2366 


29 ±1 


5.6 ±0.4 


160 ± 10 


> 450±50 


330 


1.1 


NGC 625 


29 ± 1 


1.4 ±0.1 


40. ± 2 


450±50 


720 


1.6 


NGC 784 


26 ± 1 


4.5 ±0.6 


120 ± 20 


> 450±50 


470 


1.0 


Ho II 


27 ±1 


6.5 ±0.8 


180 ± 20 


> 570±70 


820±20 


0.70±0.13 


NGC 5253 


45 ± 1 


9.0 ±0.9 


400 ± 40 


> 450±50 


450±20 


1.0±.0 


NGC 6822 


2.3 ±0.1 


3.1 ± 1.1 


7.3 ±2.5 


> 20. ±13 


140 


0.14 


NGC 4214 


43 ± 2 


3.1 ±0.9 


130 ± 40 


> 810±190 


400±10 


2.0±0.5 


NGC 1569 


21 ±2 


12 ±1 


240 ± 10 


> 450±50 


170±10 


2.6±0.3 


NGC 4449 


160 ± 10 


6.0 ±0.5 


970 ± 70 


> 450±50 


300±10 


1.5±.l 



Note. — Col. (5) Durations are lower limits for ongoing bursts. Col. (6) Timescales were calculated using ma- 
jor axis measuremen ts and maximum rotational velocity corrected for inclination from the HyperLeda database 
jPaturel et al.l 20031) with the exception of NGC 6822 whose maximum rotational velocity was obtained from 
IWeldrake et al. N2003h . Uncertainties were not reported on the rotational velocities of five galaxies, thus the dynam- 
ical timescale lacks uncertainties in these cases. 
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Table 2. Stellar Mass and Energy Calculations 





M B 


Total Stellar 




% Mass 


loff FjllfTETV of 


Galaxy 




Mass (1O 6 M0) 




Created in Burst 


Burst (erg) 


Antlia Dwarf 


-10.14 


7.3 ±0.9 


8.0 ±1.0 


3.2 ±0.5 


54.2 


UGC 9128 


-12.45 


19 ±3 


2.4 ±0.4 


24 ±5 


55.5 


UGC 4483 


-12.68 


15 ±3 


1.6 ±0.4 


26 ± 7 


55.4 


NGC 4163 


-13.75 


140 ± 40 


5.5 ±1.4 


2.8 ±0.8 


55.5 


UGC 6456 


-13.85 


68 ±27 


2.4 ±1.0 


5.1 ± 2.2 


55.3 


NGC 6789 


-14.60 


100 ± 30 


1.8 ±0.5 


3.9 ±1.4 


55.5 


NGC 4068 


-14.96 


320 ± 50 


4.1 ±0.7 


4.8 ±0.8 


56.0 


SBS 1415+437 


-15.11 


240 ± 50 


2.7 ±0.6 


20. ±4 


56.5 


DDO 165 


-15.19 


270 ± 50 


2.8 ±0.5 


19 ±4 


56.5 


IC 4662 


-15.39 


420 ± 90 


3.7 ±0.8 


5.2 ± 1.2 


56.2 


ESO 154-023 


-16.21 


650 ± 90 


2.7 ±0.4 


6.4 ±0.9 


56.5 


NGC 2366 


-16.33 


370 ± 50 


1.3 ±0.2 


12 ± 2 


56.5 


NGC 625 


-16.26 


370 ± 140 


1.5 ±0.6 


3.8 ± 1.5 


56.0 


NGC 784 


-16.78 


1100 ±100 


2.6 ±0.3 


4.7 ±0.7 


56.5 


Ho II 


-16.92 


450 ± 60 


1.0 ±0.1 


12 ± 2 


56.5 


NGC 5253 


-16.98 


2200 ± 300 


4.3 ±0.6 


4.7 ±0.7 


56.9 


NGC 6822 


-17.86 


28 ±4 




1.0 ±0.2 


53.9 


NGC 4214 


-17.02 


400 ± 100 


1.0 ±0.2 


18 ±5 


56.7 


NGC 1569 


-17.76 


1000 ± 100 


1.0 ±0.1 


5.8 ±0.7 


56.6 


NGC 4449 


-18.02 


3000 ± 500 


2.3 ±0.4 


7.0 ± 1.5 


57.2 



Note. — Col. (2) Integrated absolute magnitude of the galaxy adjusted for extinction in the 
B band. Cols. (3) The mass derived from the SFH does not incorporate all the mass from the 
optical extent of the galaxies as the field of view in some cases does not extend to the outer 
optical limits. While this technically makes the mass a lower limit, it is a small effect as most of 
the optical area is covered in each system and the amount of mass at these outer radii is small 
compared to the mass in the inner parts of the systems. The exception in NGC 6822 where the 
observations cover less than half of the optical disk. The total stellar mass is a lower limit for 
the galaxy. The mass-to-light ratio is not calculated for NGC 6822 as the mass and luminosity 
measurements cover significantly different areas of the galaxy. Col. (4) Uncertainties arc 
derived assuming errors are dominated by uncertainties in mass; uncertainties in luminosity 
are not included but are negligible at ~ 0.2 mag. Col. (5) The percent of mass created in 
each burst are likely upper limits. The high SFRs in the galaxies prevent a complete census of 
the older (t> 6 Gyr) populations. This is particularly true in the case of SBS1415±437 whose 
ancient SFH remains uncertain due to limited depth of photometry. 
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Age (Gyr) Age (Myr) 



Fie. 1. This four panel figure shows the image, photometry, and SFH results for UGC 9128. Panel (a) presents the V band 

image of UGC 9128 from the HST ACS instrument. Panel (b) shows the CMD with average photometric errors per magnitude 
bin. The MS and HeB stages of evolution are marked. The regions on the CMD labeled (A), (B), and (C) correspond to three 
different episodes of recent SF that are not spatially coincident in the galaxy. Event (A) in the CMD includes BHeB stars, 
RHcB stars, and MS stars. Events (B )and (C) include the BHeB stars highlighted in CMD. The RHeB stars from these events 
are blended with the red clump while the MS stars born simultaneously with the BHeB stars in regions (B) and (C) lie below 
the photometric depth in the CMD. The reconstructed SFH for the lifetime of UGC 9128 is presented in panel (c). The solid 
red line indicates the average SFR over the past 6 Gyr (i.e., b rece nt = 1) and the dashed red line is twice this average (i.e., 
brecent = 2). An expansion of the SFH during the last 500 Myr is presented in panel (d). The SF events (A), (B), and (C) 
identified in the CMD in panel (b) and the intervening periods of quiescent SF are reconstructed in the SFH. 
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Fig. 2. — The SFH of the last 1.5 Gyr for the twenty galaxies. The durations of the bursts 
are shaded in blue. The dynamical timescales of the galaxies are plotted in a green bar for 
comparison. The solid red lines indicate the average SFR over the last 6 Gyr for the galaxy 
(i.e., b r ecent = 1) and the dashed red line is twice this average (i.e., b recent = 2). 
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Fig. |2j — Continued 
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Fig. |2j — Continued 
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Fig. |2j — Continued 
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Fig. |2j — Continued 
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Fig. 3. — Using SF input at varying look-back times (dashed lines at 100 Myr, 200 Myr, 450 
Myr, 815 Myr), we recovered the SFHs from each single age populations (blue, red, cyan, and 
magenta lines) from photometry reaching a depth of I~ 2 mag (i.e., Antlia dwarf galaxy, top 
panel), I~ —1 mag (i.e., NGC 4068, typical of most of our sample, middle panel), and I~ —2 
mag (i.e, NGC 4214, bottom panel). Our choice in time binning at the different look-back 
times is shown as shaded gray areas. The inherent time resolution is finer at recent times 
whereas the co-variance between the time bins at intermediate and ancient times justifies a 
coarser time resolution. 
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Fig. 4. — The histogram shows the distribution of the starburst durations for the galaxy 
sample. The histogram bins are a simple scheme with each bin having the width of 100 
Myr. The actual uncertainties in the measured durations vary; the shortest durations have 
the smallest uncertainty (i.e., 450 ± 50 Myr) and the longest durations have the largest 
uncertainty (i.e., 1300±290 Myr). Also plotted is the distribution of the dynamical timescales 
for the seventeen galaxies for which rotation speeds are available. Interestingly, the duration 
of the starbursts for fifteen of the galaxies lie in the range ~ 450— ~ 650 Myr with the 
majority presenting burst duration of 450 Myr. Many of these durations are lower limits as 
the bursts are ongoing; the distribution would likely shift to even longer timescales with an 
exclusively fossil burst sample. Two galaxies in the sample show a longer duration greater 
than 1 Gyr. The timescale for the evolution of massive stars associated with "self-quenching" 
is plotted to show that the scale of these burst durations far exceed this timescale. 



